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Abstract 


The largest extinct volant birds (Pelagornis sandersi and Argentavis magnificens) and pterosaurs (Pteranodon and Quetzalcoatlus) are 
thought to have used wind-dependent soaring flight, similar to modern large birds. There are 2 types of soaring: thermal soaring, 
used by condors and frigatebirds, which involves the use of updrafts to ascend and then glide horizontally; and dynamic soaring, 
used by albatrosses, which involves the use of wind speed differences with height above the sea surface. Previous studies have sug- 
gested that P. sandersi used dynamic soaring, while A. magnificens and Quetzalcoatlus used thermal soaring. For Pteranodon, there is 
debate over whether they used dynamic or thermal soaring. However, the performance and wind speed requirements of dynamic and 
thermal soaring for these species have not yet been quantified comprehensively. We quantified these values using aerodynamic mod- 
els and compared them with that of extant birds. For dynamic soaring, we quantified maximum travel speeds and maximum upwind 
speeds. For thermal soaring, we quantified the animal’s sinking speed circling at a given radius and how far it could glide losing a 
given height. Our results confirmed those from previous studies that A. magnificens and Pteranodon used thermal soaring. Conversely, 
the results for P. sandersi and Quetzalcoatlus were contrary to those from previous studies. P. sandersi used thermal soaring, and Quet- 
zalcoatlus had a poor ability both in dynamic and thermal soaring. Our results demonstrate the need for comprehensive assessments 
of performance and required wind conditions when estimating soaring styles of extinct flying species. 
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Significance Statement: 


Quetzalcoatlus (10 m wingspan), one of the largest pterosaurs, was thought to have flown over land using updrafts like condors 
and eagles. Conversely, Pelagornis sandersi (7 m wingspan), one of the largest extinct volant birds, was thought to have flown using 
dynamic soaring like albatrosses, using differences in wind speed with height above the sea surface. In this study, we used aerody- 
namic models to comprehensively quantify soaring performances and wind requirements of these extinct species and compared 
them with extant soaring birds. We found that Quetzalcoatlus was less suited to flying in updrafts than the extant birds, and P. 
sandersi was better suited to flying in updrafts above the sea, similar to frigatebirds, rather than using albatross-like dynamic 
soaring. 


Introduction is estimated to have had a wingspan of 6-7 m (Table 1) (8, 9). The 


Flying animals have evolved a wide range of body sizes. Among 
them, there have been exceptionally large species of birds 
and pterosaurs (Fig. 1). Among the many extinct giant bird 
species (1-7), Pelagornis sandersi (from the late Oligocene, ap- 
proximately 25-28 million years ago [Ma]) and Argentavis mag- 
nificens (from the upper Miocene, approximately 6 Ma) are the 
argest volant birds. Their estimated wingspans reached 6- 
7 m (1-4), twice as large as that of the wandering albatross 
Diomedea exulans], the extant bird with the longest wingspan (Ta- 
ble 1). Several large species of pterosaurs appeared in the Cre- 
taceous period. Pteranodon, arguably the most famous pterosaur, 





azhdarchids are one of the most successful Cretaceous pterosaur 
groups and include several large species with wingspans of ap- 
proximately 9-11 m (Table 1) (10-13). Although their huge sizes 
have been led debate about whether they were flightless (14-16), 
Quetzalcoatlus northorpi, an azhdarchid species, is often regarded 
as one of the largest flying animals in history (16). 

How and how well these giant animals were able to fly has 
fascinated researchers across disciplines for over a century (23). 
This is because the question is not only interesting as a biophys- 
ical question in its own night, but also because it contributes to 
unraveling a wide range of issues such as the lifestyle of these 





OXFORD 


UNIVERSITY PRESS 





Competing Interest: The authors declare no competing interest. 

Received: September 16, 2021. Revised: February 18, 2022. Accepted: March 3, 2022 

© The Author(s) 2022. Published by Oxford University Press on behalf of the National Academy of Sciences. This is an Open Access article 
distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/), which permits 
unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited. 


zzoz Aew 9z uo 1senBb Aq ,oz9rs9/Eezooe6d/ 1/1 /ejonse/snxeuseud/Woo'dno'olwepeoe//:sdiy Wo pepeojumoq 


2 | PNAS Nexus, 2022, Vol. 1, No. 1 





California condor 


sal 





wandering albatross 


mg 


magnificent frigatebird 






Pteranodon 


kori bustard 





Argentavis magnificens 





Pelagornis sandersi 


Quetzalcoatlus r— 
1m 


Thermal soarer Leas Dynamic soarer e Poor soarer 


Fig. 1. A size comparison and soaring styles of extinct giant birds (P. sandersi and A. magnificens), pterosaurs (Pteranodon and Quetzalcoatlus), the largest 
extant dynamic soaring bird (wandering albatross), the largest extant thermal soaring terrestrial bird (California condor), a large extant thermal 
soaring seabird (magnificent frigatebird), and the heaviest extant volant bird (kori bustard). The icons indicate dynamic soarer, thermal soarer, and 
poor soarer, and summarize the main results of this study. The pink arrows indicate the transition from a previous expectation or hypothesis to the 


knowledge updated in this study. 


species, their role in paleoecosystems, and the drivers of morpho- 
logical evolution, diversification, and extinction of giant species 
over geological time (6, 16, 24-26). Their huge size must have sig- 
nificantly affected their flight because, with increasing size, the 
power required to fly increases faster than the power muscles can 


produce via the flapping of wings (27-29). Hence, this physical con- 
straint has resulted in 2 debated arguments about the flight of 
extinct giants. The first debate is about whether and how they 
were able to take off (4, 14, 30, 31). The second debate concerns 
their soaring flight style. Due to the high costs of flapping that 
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Table 1 Morphological values of examined species. 




















Wing area Wing loading 
Species Mass [kg] Wingspan [m] [m?] Aspect ratio [N/m?] Ref 
; ‘ ‘ . 6.06, 6.13, 6.40, 13.0, 14.0, and 51.0-87.4* and 
Extinct animals Pelagornis sandersi 21.8 and 40.1 and 7.38 2.45-4.19 5.0 93.9-1610 (1) 
Argentavis magnificens 70.0 7.00 8.11 6.04 84.7 (4) 
ee 36.7 5.96 1.99 7.9 181 (9) 
18.6 5.34 2.26 2.6 80.7 (15) 
Quetzalcoatlus 259 9.64 11.4 8.18 224 (9) 
Dynamic soaring wandering albatross 8.64 3.05 0.606 5.4 140 (17)° 
birds black-browed albatross 3.55 2.25 0.376 3.4 87.5 (18)° 
white-chinned petrel 1.37 1.40 0.169 1.6 795 (19) 
Thermal soaring magnificent frigatebird 1/52 2.29 0.408 2.8 36.5 (20) 
birds California condor 9.50 2.74 1.32 5.70 70.6 (4) 
brown pelican 2.65 2.10 0.450 9.80 57.8 (20) 
black vulture 1.82 1.38 0.327 5.82 54.6 (20) 
white stork 3.40 2.18 0.540 7.42 61.8 (4) 
Non-soaring bird kori bustard 11.9 2.47 1.06 5.76 110 (21)4 
Motor glider Schleicher ASK 14 340 14.3 12.6 16.2 265 (22) 





®With 21.8 kg mass. 
>With 40.1 kg mass. 


°We used the averages calculated from the morphological values of males and females in the cited references. 


4We used the morphology data available in the Flight program. 


stems from their large body size, large extant birds prefer to fly 
utilizing wind energy or convection, that is, they prefer to soar 
(14, 32). Hence, it is presumed that extinct large animals also em- 
ployed soaring flight as their primary mode of transportation (1, 
4, 16). However, the kind of soaring flight style they employed is 
debated (1, 4, 16, 33, 34). The present study focuses on this second 
debate. 

There are 2 main soaring flight styles among extant birds: dy- 
namic soaring and thermal (and slope) soaring (21). In dynamic 
soaring, birds extract flight energy from wind shear—the verti- 
cal gradient in horizontal wind speed over the ocean (Fig. 2A- 
B). Extant seabirds (e.g. albatrosses, shearwaters, and petrels) 
employ this soaring style and can routinely travel hundreds of 
kilometers per day over the sea. In thermal soaring, birds first 
fly circling in warm rising-air columns (thermals). They climb 
to a substantial height and then glide off in the desired direction 
while losing their height (Fig. 2C-E). By repeating this up-down 
process, birds travel over vast distances. Various terrestrial bird 
species (e.g. vultures, eagles, and storks) and seabirds (e.g. frigate- 
birds and pelicans) employ thermal soaring (20). Terrestrial birds 
that utilize thermal soaring are also able to use the upward wind 
deflected by a mountain or hillside instead of the thermal for their 
ascent (32); this is termed “slope soaring.” Another form of slope 
soaring called wave-slope soaring uses the upward wind gener- 
ated by waves (35). Wave-slope soaring is thought to be used 
by both thermal soaring seabirds (pelicans) and dynamic soar- 
ing seabirds (albatrosses) (35), but since no bird exclusively uses 
wave-slope soaring, we did not include it in the main soaring 
styles here. 

Previous studies estimated that P. sandersi was a dynamic soarer 
(1), and A. magnificens and Quetzalcoatlus were thermal soarers (4, 
16). Pteranodon was estimated to be a dynamic soarer (16), but 
subsequent studies have indicated that it was a thermal soarer 
(34, 36). See Quantification of soaring styles in previous studies in Ma- 
terials and Methods for details about previous studies on this 
topic. 








To estimate the potential soaring styles of these extinct ani- 

mals, it is necessary to quantify their soaring performance, e.g. 
potential speed and efficiency of soaring, as well as the required 
wind speed to sustain soaring flight. Valuable indicators of dy- 
namic soaring performance are the maximum travel speed and 
the maximum upwind speed averaged over 1 dynamic soaring 
cycle (Fig. 2B) (37-39). Additionally, it is necessary to evaluate 
the minimum horizontal wind speed required for sustainable dy- 
namic soaring (40). Thermal (and slope) soaring performance is 
well-quantified by 2 indicators: the glide ratio, i.e. the ratio of 
the distance the animal traverses to the height the bird loses to 
cover that distance in the gliding phase (Fig. 2E), and the sink- 
ing speed of the animal circling in a given radius during the up- 
ward soaring phase (Fig. 2D). This sinking speed during circling 
corresponds to the upward wind speed required to ascend in a 
thermal. Because thermals have a stronger updraft in the cen- 
ter (Fig. 2C), the animal needs to achieve not only low sinking 
speed but also a narrow circle radius to efficiently ascend using 
a thermal. For slope soaring, the minimum value of sinking speed 
(sinking speed at infinite radius) would be an appropriate met- 
c for performance, as the circle radius is less restrictive in slope 
oaring. 
However, the soaring performances and required wind condi- 
tions have not been comprehensively evaluated for P. sandersi, 
A. magnificens, Pteranodon, and Quetzalcoatlus (summarized in Ta- 
ble 2). A total of 3 knowledge gaps are highlighted in Table 2: 
(1) the dynamic soaring performance and the minimum required 
wind speed have never been evaluated; (2) the thermal soaring 
performance in the soaring up phase and the required updraft 
wind speed have not been evaluated for P. sandersi, Pteranodon, 
and Quetzalcoatlus; and (3) despite recent studies showing that 
the body masses of Pteranodon and Quetzalcoatlus were approxi- 
mately 3 times heavier than previously expected (9, 14, 15), and 
that pterosaurs’ wings had a higher profile drag than that of birds 
(34), the soaring performances of these new heavy body masses 
and higher drags have rarely been evaluated. 








oO 4, 
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Fig. 2. Schematics of dynamic soaring and thermal soaring. (A) Example of a 3D track of dynamic soaring. Dynamic soaring species repeat an up and 
down process with a shallow S-shaped trajectory at the sea surface. By utilizing wind gradients, a species can fly without flapping. (B) Example of a 2D 
dynamic soaring trajectory of 1 soaring cycle. The travel speed averaged over 1 cycle is defined as the travel distance in 1 cycle (d) divided by the 
soaring period, and the upwind speed averaged over 1 cycle is defined as the upwind travel distance in 1 cycle (dup) divided by the soaring period. (C) 
Schematic of a thermal soaring cycle. (D) In the soaring up phase, a species soars in a steady circle. When there is upward wind that is greater than a 
species’ sinking speed, the species can ascend in the thermal. The upward wind is stronger in the center of a thermal; therefore, achieving a small 
circle radius is advantageous for thermal soaring. (E) In the gliding phase, a species glides in a straight line. The rate of horizontal speed to the sinking 


speed is equal to the rate of horizontal distance traveled to the height lost. 


In this study, we aimed to address these knowledge gaps and 
identify the potential soaring styles of these extinct giant birds 
and pterosaurs. To this end, we used physical models and recent 
morphology estimates to quantify the performance and wind con- 
ditions required for dynamic and thermal soaring in these extinct 


animals, and compared them with those of extant soaring birds. 
See Models in Materials and Methods for details about the em- 
ployed models and parameter values. 

Another aim of this study is to provide a framework for the ob- 
jective evaluation of the soaring performance of extinct animals 
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Table 2 Summary of previous studies. Previous studies that quantified the soaring performances and required wind conditions of P. 
sandersi, A. magnificens, Pteranodon, and Quetzalcoatlus with recent heavy body mass estimates. 
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Pelagornis sandersi Dynamic soaring (1) _— 


Glide polar (1) oa 


Glide polar (1) — 





Argentavis magnificens Thermal soaring (4) _ —_— 


Circling envelope (4) Glide polar (4) Circling envelope (4) 





Thermal soaring (34) 
Pteranodon : : — _— 
Dynamic soaring (16)! 


Quetzalcoatlus Thermal soaring (16)! _ — 


= Glide polar (34) — 





TA principal component analysis (PCA) using 3 morphological information (logarithm of each of the animal’s weight, wing area, and wingspan) has been performed 
in the previous studies for birds (41), bats (42), and pterosaurs (9), respectively. In the previous study (16), these second and third principal components were 
compared among birds, bats, and pterosaurs to estimate the soaring styles of pterosaurs. 





using models based on physical principles that all animals adhere 
to and serve as a stepping stone for invoking further research on 
this topic. Therefore, the assumptions and simplifications in our 
model, as well as our results and their implications, are impor- 
tant outcomes for future research and are explored in depth in 
the Discussion section. 








Results 
Dynamic soaring 


We quantified the dynamic soaring performance (Fig. 2B) and re- 
quired wind speeds using a physical model and a numerical op- 
timization method. This method has been developed in the en- 
gineering field and provides a framework to quantify dynamic 
soaring performances and required wind conditions for gliders 
and birds (37, 40, 43, 44). However, despite its effectiveness, the 
only animal to which this technique has been applied is the wan- 
dering albatross (40, 44); it has never been applied to extinct giant 
flyers. We applied this framework to the 4 giant extinct species 
and 3 extant dynamic soaring bird species with various sizes rang- 
ing from 1 to 9 kg (i.e. the white-chinned petrel [Procellaria aequinoc- 
tialis], the black-browed albatross [Thalassarche melanophris], and 
the wandering albatross). As the exact shape of the wind gradient 
remains poorly understood, we conducted the calculation under 7 
different wind conditions (Fig. 3A—-C) (40, 44). In addition, we added 
an important modification to the previous models: the animal’s 
wings do not touch the sea surface during their flight (Fig. 3D; 
and see Eq. 17 and its description in Materials and Methods for 
details). 

Our computation results indicate that A. magnificens, Pteran- 
odon, and Quetzalcoatlus showed lower dynamic soaring perfor- 
mances and higher required wind speeds for dynamic soaring 
than the extant dynamic soaring species under all wind condi- 
tions tested in this study (Fig. 4; Figures $1 and S2, Supplemen- 
tary Material). Therefore, these species could not have employed 
dynamic soaring. 

For P. sandersi, the results were highly dependent on body mass. 
With heavy body mass estimates (40.1 kg), P. sandersi required 
higher wind speeds than extant dynamic soaring species, irre- 
spective of the wind conditions. The performance was superior 
to extant species for some morphology estimates when the shear 
height was far from the sea surface (sigmoidal wind condition 
with hy = 5m), but inferior when the wind speed change was 
located close to the sea surface (logarithmic model and sigmoidal 
wind condition with hy = 1 and 3 m). In particular, their up- 
wind speed was inferior to that of extant species, resulting in an 
inefficient flight at their upwind destination and, consequently, 














Log Sigmoid 
A B Cc 
1 ae 5 = 7/6 6 = 3/6 5 
£38 h h 
ce y 
2 4 
Io5 z 





Fig. 3. Wind shear models explored in this study. (A) Logarithmic wind 
gradient model. The wind speed at height 10 m was defined as Wp. (B) 
Sigmoidal wind shear model with a wind shear thickness of 7 m (5 = 7/6) 
and a shear height (hw) of 1, 3, or 5 m. (C) Sigmoidal wind shear model 
with a wind shear thickness of 3 m (5 = 3/6) and a shear height of 1, 3, or 
5m. The maximum wind speed of the sigmoidal model is represented as 
Wmax- (D) Schematic of a soaring bird. Its height from the sea surface is 
represented as z and the height of the wingtip is represented as Zying. 
We constrained the models so that the wing tip did not touch the sea 
surface, 1.€. Zwing > 0. 


limiting their travel direction. When lower body mass estimates 
were used (21.8 kg), P. sandersi required lower wind speeds, but its 
performance was distinctively lower than that of extant species. 
Hence, P. sandersi required harsh wind conditions for dynamic 
soaring when a 40.1 kg body mass was assumed, and it was poor 
at dynamic soaring when a 21.8 kg body mass was assumed. 

The performances of all species varied with the value of hy, and 
the variation was especially distinct for large species in contrast 
to that of white-chinned petrels (Fig. 4; Figures S1 and S2, Sup- 
plementary Material). This variation was due to wingtip bound- 
ary conditions (Fig. 3D and Eq. 17). Animals can attain more en- 
ergy when passing through large wind speed gradients, but when 
large gradient changes are close to the sea level, large animals are 
unable to use the wind speed gradient efficiently because their 
wings limit the altitude available to them. Although the long, thin 
wings that reduce drag in extant dynamic soaring birds are suited 
for dynamic soaring (21, 45), our detailed dynamic models have 
shown that excessively long wings can also inhibit efficient dy- 
namic soaring. 


Thermal soaring 


The thermal soaring performances and the required upward wind 
speeds (Fig. 2D and E) were quantified using the established 
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Fig. 4. Minimum required wind speeds and dynamic soaring performances of extinct and extant animals. (A) Results of the logarithmic wind model. 
(B) Results of the sigmoidal wind model with a wind shear thickness of 7 m (5 = 7/6) and wind shear height (hw) of 3 m. The first row shows the 
minimum required wind speed for sustainable dynamic soaring. The second row shows the maximum travel speed averaged over 1 soaring cycle, in 
response to wind speed. The third row shows the maximum upwind speed averaged over 1 soaring cycle, in response to wind speed. 


framework, i.e. glide polars and circling envelopes. A glide po- 
lar is a graphical plot of the sinking speed vs. horizontal speed 
when a bird glides in a straight line (21). We can determine the 
maximum glide ratio and the associated travel speed of flyers by 
identifying the line that passes through the origin and tangents of 
the glide polar plot. The inverse of the line slope and the speed at 
the tangent point correspond to the maximum glide ratio and the 
associated horizontal speed, respectively. The circling envelope is 
a graphical plot of the sinking speed (1.e. equivalent to the required 
upward wind speed for ascent in a thermal) vs. the radius of a turn 
when a bird glides in a steady circle (21, 20). 

We quantified the thermal soaring performances and the re- 
quired upward wind speeds for the 4 extinct species, 5 ex- 
tant thermal soaring species (the magnificent frigatebird [Fregata 
magnificens], the black vulture [Coragyps atratus], the brown peli- 
can [Pelecanus occidentalis], the white stork [Ciconia ciconia], and the 
California condor [Gymnogyps californianus]), and the kori bustard 
(Ardeotis kori), the heaviest extant volant bird that does not soar. 
The performances reported in a previous study of the Schleicher 
ASK-14, a motor glider with a 14 m wingspan similar to that of 
Quetzalcoatlus, are also presented for comparison (22). 

All extinct species showed high gliding performances with 
maximum glide ratios ranging from 10 to 17 (Fig. 5A), which are 
comparative to those of extant species. 

With respect to the soaring up phase, all of the extinct giant 
flyers, except for Quetzalcoatlus, had performances equivalent to 








or better than the extant soaring species (Fig. 5B). As shown pre- 
viously (4), the circling up performance of A. magnificens was com- 
parable to that of the California condor, one of the largest living 
thermal soarers. The performance of Pteranodon was comparable 
to that of living thermal soarers; however, it was rather low when 
the higher body mass estimate (36.7 kg) was employed. The ther- 
mal soaring ability of P. sandersi when a light mass was assumed 
21.9 kg) was outstanding. It outperformed several extant thermal 
soaring species in soaring up ability, and was even comparable 
to the magnificent frigatebirds, the champion of thermal soaring 
among extant species. Even with a heavier body mass estimate 
40.1 kg), P. sandersi still outperformed or was comparable to sev- 
eral other species. 

Among the 4 extinct giant animals investigated in this study, 
the soaring up performance of Quetzalcoatlus was exceptionally 
ow. It required the strongest upward wind speed and the widest 
circle radius. Its performance was even lower than that of the 
kori bustard, one of the heaviest volant extant bird species, which 
spends most of its time on land and only fly in emergencies, such 
as when under predation risk. 

Figure 5 also demonstrates how gliders and animals differ with 
regards to what they prioritize in thermal soaring. A low sink- 
ing speed and a narrow circle radius are advantageous during 
the ascent phase of thermal soaring, whereas a high glide ra- 
tio and flight speed are advantageous during the glide phase, 
the benefits of which are not simultaneously satisfied. Natural 
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Fig. 5. Glide polars (A) and circling envelopes (B) of extinct species, extant thermal soaring species, and the kori bustard, the heaviest rarely flying bird. 
In (A), the maximum glide ratios of each species are shown on the right side of species names. Points represent the horizontal speed and sinking speed 
at the maximum glide ratio of each species. Gray lines represent glide ratios (5, 10, 15, 20, 25, and 30). (B) Shows circling envelopes defined as the 
minimum sinking speed against the circle radius. The left end of the curve is for a bank angle of 40°. The smaller the bank angle, the larger the circle 
radius. A linear wingspan reduction is assumed for birds and a fixed wingspan is assumed for pterosaurs. The lift coefficient of circling envelope (C7) is 
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indicates that while its high flight speed results in a larger cir- 
cle radius than that of animals, its maximum glide ratio is nearly 
30, far superior to that of birds and pterosaurs. Some modern 
sailplanes (such as Eta) have glide ratios of up to 70 and aspect ra- 
tios up to 50 (47). Therefore, the distinct difference in performance 
between the glider and animals (Fig. 5) indicates the difference in 
“selection pressure” and constraints on morphology (i.e. materials 
available for the wing) between them. 





Discussion 


Although several previous studies have investigated the soaring 
performance of extinct species, there have been several evalua- 
tion gaps. In the present study, we filled these gaps using physical 
models of soaring birds. We computed and compared the dynamic 
and thermal soaring performances and the required wind condi- 
tions for soaring of 4 extinct giant flyers with those of extant dy- 
namic and thermal soaring species, which enabled us to examine 
the soaring style of extinct giants from multiple perspectives. Our 
results indicate that A. magnificens was a thermal soarer, confirm- 
ing a previous study (4). For Pteranodon, there were 2 previous com- 
peting claims: dynamic soarer (16) and thermal soarer (34). Our re- 
sults supported that they were thermal soarers. For Quetzalcoatlus 
and P. sandersi, our results were in contrast to previous studies. Our 
results indicate that Quetzalcoatlus could not efficiently perform 
dynamic nor thermal soaring. In addition, although P. sandersi was 
considered a dynamic soaring species in a previous study (1), our 
results suggest that it was a thermal soaring bird. We discuss our 
results in detail for Pteranodon, Quetzalcoatlus, and P. sandersi, and 
then describe future issues that need to be addressed for a better 
understanding of the soaring styles of extinct giant species. 





Pteranodon 


There have been conflicting estimates of the soaring style of Pter- 
anodon. It was claimed that they used thermal soaring based on 
quantitative evaluations of glide polars obtained from wind tun- 
nel experiments using reconstructed wings of pterosaurs (34) and 
of the strength of the pterosaur’s wing membranes using a physi- 
cal model (36). In contrast, it was claimed that they used dynamic 
soaring based on a comparison of morphology between birds, bats, 
and pterosaurs with a principal component analysis (PCA), which 
showed that Pteranodon and albatrosses were located in similar re- 
gions of the principal components space (16). 

However, it should be noted that the comparative PCA method 
employed in the latter study does not provide quantitative soar- 
ing performances, as cautioned by the proposer of the original 
method (41) (this issue will be discussed in detail in the next sec- 
tion, “Quetzalcoatlus”). Hence, there has been no direct evaluation 
of the dynamic soaring performance of Pteranodon. Our results 
fill this knowledge gap and show that Pteranodon was not good at 
dynamic soaring and had a thermal soaring ability comparable to 
that of extant species, as has already been shown (34). 

Based on these points, and as many Pteranodon have been found 
from deposits in the center of the Western Interior Seaway (48), we 
suggest that they travelled over oceans by thermal soaring like 
frigatebirds. 








Quetzalcoatlus 

There has been a heated debate about the flight capability of Quet- 
zalcoatlus. The focal issue has been whether or not Quetzalcoatlus 
could take off. Researchers are divided between the opinion that 
it was too heavy to take off (14,15) and the opinion that it was 


able to take off by using quadrupedal launching, like some bats 
(16, 30) (but also see the recent discussion on the launching of 
Quetzalcoatlus (49)). In addition, detailed observations of fossils 
are also presented as evidence that the giant azhdarchids, in- 
cluding Quetzalcoatlus, were capable of flight; for example, a huge 
deltopectoral crest on their humeri, which would have anchored 
muscles for flapping flight (50). 

Although there is some debate as to whether or not giant 

pterosaurs could have taken off, it has been widely accepted that 
if they were able to take off their primary mode of travel would 
have been thermal soaring rather than flapping flight (16). Wit- 
ton and Habib reported that the flapping flight of Quetzalcoatlus 
required anaerobic movement and was difficult to sustain for a 
ong period; therefore, it must have relied on wind energy for long- 
distance travel (16). Based on a comparison of morphology be- 
tween birds, bats, and pterosaurs with a PCA, these authors also 
concluded that Quetzalcoatlus used thermal soaring (9, 16). 
Our results revealed that Quetzalcoatlus had a poor performance 
to ascend using thermal or deflected upward wind. It required a 
arger circle radius and stronger updraft than the extant thermal 
soaring species, and even than the kori bustard that does not soar 
and spent most of their time on land. This suggests that the wind 
conditions under which Quetzalcoatlus could conduct sustainable 
thermal soaring were limited. The ability to climb in weak up- 
drafts at low altitudes is particularly important for sustainable 
thermal soaring. At high altitudes, strong thermals can create up- 
drafts of 5 m/s or more, and thus the sink rate of approximately 
2 m/s of Quetzalcoatlus would have been sufficient to exploit such 
strong updrafts. However, it should be noted that thermals are 
generally weaker and narrower at low altitudes (51). Therefore, ev- 
ery time an animal starts to climb a thermal, they are obliged to 
soar in the poor thermal conditions at low altitudes. For example, 
the largest existing thermal soaring species can only climb at low 
speeds at low altitudes because of the weaker updrafts of ther- 
mals and the bird’s high sink rate involved in reducing the circle 
radius to accommodate the narrower thermal width (51). Hence, 
the poorer thermal soaring ability of Quetzalcoatlus than that of ex- 
tant species would have greatly reduced the number of thermals 
they could have climbed, which would have made it difficult for 
them to stay aloft long enough to find a strong thermal. 

The poor thermal and slope soaring performance of Quetzal- 
coatlus was due to the large wing loading associated with the large 
body size. As shown in Materials and Methods, the circle radius 
was proportional to the wing loading (Eq. 21), and the minimum 
sinking speed during straight glide was proportional to the wing 
loading to the power of one-half, when the effect of the wing 
length adjustment and the lift coefficient dependence of the pro- 
file drag coefficient were negligible (Eq. 20). Since the wing loading 
increases with body mass, a giant Quetzalcoatlus required thermals 
with a wider radius and stronger updraft for thermal and slope 
soaring. 

The wing loading also explains why the results of the present 
study are not consistent with the claims of previous studies 
that Quetzalcoatlus was adapted to thermal soaring (9, 16). In the 
previous study, soaring ability was assessed from 2 variables re- 
lated to wing loading and aspect ratio with PCA (second and 
third principal components obtained by PCA on the logarithms 
of body weight, wing area, and wingspan; see (41, 42) for details 
of the method and data). Note that these variables are not ex- 
actly equal to wing loading and aspect ratio. In particular, the size- 
dependence has been removed from the principal component 
related to wing loading. However, thermal soaring performance 
is inevitably size-dependent and, therefore, using performance 
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and wind requirements calculated from morphology based on the 
aws of physics (as conducted in this study and (26, 45)) are more 
accurate. 

Interestingly, a recent research on Quetzalcoatlus suggested that 
their wing loading is even higher than previously estimated (49). 
The authors suggested narrower pterodactyloid wings, which are 
ittle or not attached to the hind limbs, while many previous stud- 
ies adopted the wings that are attached to the hind limbs and have 
a broad, bat-like shape (16, 36, 52). Although no specific morpho- 
ogical estimates are provided in the study (49), if their argument 
is correct, Quetzalcoatlus would have had a smaller wing area, and 
thus even greater wing loading than previously estimated, result- 
ing in less thermal soaring capability than that estimated in our 
study. 

Anatomical studies of the azhdarchid pterosaurs have reported 
that their skeletal structure shows adaptations to terrestrial walk- 
ing and suggested that they were terrestrial foragers (52, 53). Fur- 
thermore, a recent phylogenetic analysis showed that the azh- 
darchoid pterosaurs differed from other pterosaurs, in that they 
had evolved in a manner that increased the cost of transport for 
flapping flight and the sinking speed of gliding (26). Taking into 
account the adaptations for walking (52, 53), the humeri feature 
indicating flapping flight capability (50) but not sustainable flap- 
ping flight (16), the phylogenetic tendency of decreasing flight ef- 
ficiency (26), and the low thermal soaring ability shown here, we 
suggest that the flight styles of Quetzalcoatlus and other similar- 
sized azhdarchid species were similar to those of the bustard or 
ground hornbill that are short-range flyer and spend most of their 
time on land. Several previous studies overstated the flight abil- 
ity of large azhdarchid pterosaurs such that they had a “capacity 
for nonstop flights exceeding 10,000 miles” (54) and an “unprece- 
dented capacity for trans-continental flight" (Fig. 8 of (55)). How- 
ever, our results indicate that these claims seem highly unlikely. 











Pelagornis sandersi 


As P. sandersi was found close to the coast, this species is thought 
to have lived an oceanic existence by soaring over the sea (1). Pre- 
viously, it has been reported that P. sandersi was a dynamic soarer 
ike the albatross, rather than a thermal soarer like frigatebirds 
1). However, we argue that this species is highly adapted for ther- 
mal soaring. The conclusion of the previous study was based on 
the glide polars of P. sandersi, which were more similar to those of 
the wandering albatrosses than those of frigatebirds; glide perfor- 
mance was the only criterion used to evaluate its soaring style. 
n this study, we quantified other performances and the required 
wind conditions, which enabled us to evaluate the soaring style of 
P. sandersi from multiple perspectives. 

Our results indicated that the dynamic soaring performance 
of P. sandersi was generally inferior to that of extant dynamic 
soaring species, although there were substantial variations de- 
pending on wind conditions and morphology estimates. One of 
the factors contributing to the poor dynamic soaring ability of 
this species was an inability to efficiently exploit the wind speed 
gradient due to long wings limiting the height above sea level at 
which the bird could fly. This effect could not be assessed using 
the glide polars alone. 

Conversely, the thermal soaring ability of P. sandersi was out- 
standing, regardless of the morphology estimates used. The ther- 
mal soaring ability of an animal is largely dependent on its wing 
loading, as previously discussed. Therefore, the reason why P. 
sandersi showed such a high performance is because of its low 
wing loading despite its huge size. Considering that P. sandersi was 











found close to the coast, this species is expected to have been 
well-adapted to capture weak updrafts above the sea by using 
thermal soaring, and that it was able to stay aloft for a long period 
of time with limited flapping and traveled long distances, similar 
to frigatebirds (56). 


Future issues 


In this section, we discuss some of the simplifications used in this 
study and issues that we believe need to be addressed in the fu- 
ture. Below, we explore 3 issues with modeling flight: flight stabil- 
ity, wind environments, and discrepancy between model predic- 
tions and actual dynamic soaring speeds. 

The first issue is flight stability. In this study, a steady wind en- 

vironment was assumed, but actual wind environments fluctu- 
ate. In such a fluctuating real-world environment, stability is an 
important factor that determines the success or failure of flight 
(57-59). To simplify our calculations, we did not address stability, 
but it is important to examine the flight stability of these extinct 
and extant birds using more detailed morphological information 
in the future. 
The second issue is that the actual wind environment experi- 
enced by both extinct and extant animals is still largely unknown. 
For dynamic soaring, the specific form of the wind speed gradi- 
ent experienced by birds is unknown—for example, whether there 
is a logarithmic or sigmoidal gust in the shadows of waves (44, 60). 
For this reason, we evaluated performance under various wind 
conditions (Fig. 3). For thermal soaring, it is also unknown how 
much updraft animals experience at a given circle radius or the 
distance between thermals. Recent advances in tracking technol- 
ogy have made it possible to record details of the motion of birds in 
dynamic and thermal soaring (51, 61-64). These data will provide 
information of the real wind environment experienced by soaring 
animals (64-68). 
tis also important to consider the palecenvironmental aspects 
of the wind environment at the time of the extinct species’ inhab- 
itation. For example, we showed that Quetzalcoatlus had a lower 
thermal soaring capacity than the extinct species. More detailed 
paleoclimatic estimates may help us to understand whether the 
species experienced a quite extreme wind environment to en- 
able them to use thermal soaring as its primary mode of trans- 
port, even with their poor soaring up ability. As another ex- 
ample, we showed that Pteranodon and Pelagornis sandershi were 
oceanic thermal soarers like frigatebirds. Frigatebirds are dis- 
tributed in trade wind zones with abundant updrafts (20, 56); 
however, Pteranodon and Pelagornis had a wider distribution (1, 53). 
These mismatches between extant and extinct oceanic thermal 
soaring species could reflect the difference in the climatic en- 
vironment between the modern era and the era of the extinct 
giants. 

Third, it should be noted that in our model of dynamic soar- 
ing, the maximum travel speed the animal could achieve was very 
high. Depending on the shape of the wind speed gradient assumed 
in our model, the animals reached maximum travel speeds of 
over 100 km/h even at realistic wind speeds of 5-10 m/s (Fig. 4), 
however, the average speeds reported using GPS for albatrosses 
and shearwaters were about 30-55 km/h (67, 69) (although, at 
wind speeds of around 20 m/s, an albatross was reported to travel 
at speeds of over 110 km/h for 9 hours (70)). There are 2 potential 
reasons for the discrepancy between this model and reality. 

The first reason is that the actual shape of the wind speed gra- 
dient is likely to be more gradual than that assumed in the present 
study. It can be seen from our results in a sigmoidal wind gradient 
(Figures S1 and S2, Supplementary Material) that the gentler the 
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wind speed variation (i.e. the larger the 5), the slower the maxi- 
mum travel speed. 

The second reason is that the upper limit of the bank angle 
change speed has not been considered. In dynamic soaring, there 
is a risk of wing failure due to the dynamic load on wings, which 
is higher than the static load during stable gliding as in thermal 
soaring. The risk of wing failure can be quantified by the ratio 
of the bending moment on the wing bones induced by lift to the 
stress that the bones can withstand. This ratio is called the “fac- 
tor of safety” or “safety factor” (71-73). Animals should avoid sit- 
uations with a low safety factor, limiting their bank angle change 
speed, and hence their flight speed. Few studies have quantified 
the safety factors of flying animals (72, 74), especially for extinct 
species, in part due to the need for detailed 3D geometrical data of 
bones and wings. However, Palmer’s doctoral thesis (75) provides a 
plausible approach. Palmer quantified the safety factor along the 
wingspan of Pteranodon by applying the beam theory and vortex 
attice method to the Pteranodon wing spar structural model. The 
computed safety factor was 2.5 at the first wing phalanx when 
Witton’s mass estimates (9) were employed. Palmer concluded 
that this value is very low for a natural structure subjected to dy- 
namic loading (73), suggesting that the dynamic load demand of 
dynamic soaring cannot be sustained by a Pteranodon wing. Apply- 
ing Palmer’s framework to other species would make for intriguing 
future research. 

In view of the above 2 points, the dynamic soaring performance 
and the required wind conditions obtained from the present nu- 
merical calculation should not be taken literally by the values 
themselves. However, since all species were evaluated under the 
same assumptions, the relative values are meaningful indicators 
for the purpose of estimating soaring style by interspecies com- 
parison, as was done in this study. A more refined analysis incor- 
porating constraints on body rolling speed will be an interesting 
challenge in the future. For this purpose, actual measurements of 
bank angles in dynamic soaring birds and assessment of the wing 
bone structure and the lift force distributions along the wingspan 
in extant and extinct species will provide important information. 

Understanding the soaring performance of extinct animals 
is an interdisciplinary issue. It is, therefore, essential to have a 
place for objective discussion between researchers with different 
backgrounds to bring their knowledge together, but such a plat- 
form seems to have been lacking in the past. In our approach, 
i.e. the physical model and the framework for comprehensively 
evaluating soaring performance, we have tried to clarify what 
assumptions were needed regarding the mechanics and the 
morphology of the animals. On the basis of our theoretical 
framework, it should, therefore, be easy for experts from var- 
ious disciplines, including paleontologists, paleoclimatologists, 
engineers, and ornithologists, to examine the validity of the 
assumptions, examine new information, such as updated mor- 
phological data, and improve the model. We hope that our theo- 
retical framework presented in this study will inspire researchers 
from various disciplines to work together to understand the soar- 
ing performance of extinct animals. 















































Materials and Methods 

Quantification of soaring styles in previous 
studies 

This section reviews previous studies on the soaring performance 
of extinct giant animals. In particular, we focus on which indices 
were quantitatively evaluated for each species. 


Pelagornis sandersi 


P. sandersi is predicted to be a dynamic soarer rather than a ther- 
mal soarer as its glide polar (and glide ratio that can be derived 
from its glide polar) is more similar to those of living dynamic 
soarers than those of living thermal soarers (1). However, this 
means the understanding of this species’ soaring style has been 
based on just 1 metric, that is, its glide ratio (Table 2). Hence, eval- 
uating other metrics of this species could provide a more accu- 
rate estimate of the soaring style of this species. A previous study 
cautiously calculated the glide polars of P. sandersi for 24 combi- 
nations of estimates (Table 1) to deal with morphological uncer- 
tainty. Hence, we also employed these estimates in this study. 








Argentavis magnificens 


A. magnificens is expected to be a thermal soarer. A previous study 
reported that the thermal soaring performance and required wind 
conditions of this species were comparable to living thermal soar- 
ing species based on glide polars and circling envelopes (4). This 
result is consistent with the fact that an Argentavis specimen was 
found on the foothills and Pampas of Argentina, far from coast- 
lines (4). 


Pteranodon and Quetzalcoatlus 


Although assessments of the soaring abilities of Pteranodon and 
Quetzalcoatlus have been a long-standing issue, there is still lack 
of a comprehensive understanding of their soaring style due to 
several uncertainties in the estimates of their morphology, es- 
pecially because of the significant changes in weight estimates 
around 2010. Previously, it was estimated that Pteranodon had a 
wingspan of around 7 m and a body mass of 16 kg, while Quet- 
zalcoatlus had a wingspan of around 11 m and a body mass of 
50-70 kg. Based on these estimates, previous studies argued that 
they were adapted to thermal soaring (76, 77) and others argued 
that they could also employ dynamic soaring (33). Around 2010, 
however, several studies with different approaches suggested that 
pterosaurs were much heavier than previously expected (9, 14, 
15) (Table 1). For example, Witton estimated that Pteranodon was 
36.7 kg with a 6.0 m wingspan, and Quetzalcoatlus was 259 kg with 
a 9.6 m wingspan (9). Henderson estimated that Pteranodon was 
18.6 kg with a 5.3 m wingspan, and Quetzalcoatlus was 544 kg with 
a 11.2 m wingspan (15). Witton and Habib argued that Pteranodon 
was a dynamic soarer and Quetzalcoatlus was a thermal soarer by 
comparing their morphology with those of extant soaring animals 
using PCA (9, 16, 41). Conversely, a recent study quantified the 
cost of transport and sinking speeds during gliding in pterosaur 
species and showed that azhdarchoid pterosaurs, including Quet- 
zalcoatlus, had lower flight efficiency than the other pterosaurs 
(26). Despite these studies, the performance and wind require- 
ments of dynamic and thermal soaring in these species have not 
been comprehensively quantified. 

Furthermore, pterosaurs have a wing morphology that is com- 
pletely different from that of birds and bats. Some studies re- 
ported that the wings of pterosaurs would have been associated 
with high-profile drag (34, 78). Palmer experimentally measured 
profile drag in a wind tunnel experiment using reconstructed 
pterosaur wings, and determined glide polars of Pteranodon-sized 
pterosaurs with various body mass estimates including the most 
recent heavy estimates (34). Palmer concluded that Pteranodon had 
the comparative thermal soaring ability to that of extant therma 
soaring species, and that they adopted a slow flight speed. This 
conclusion was further reinforced by a subsequent work, which 
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quantitatively showed that the pterosaur’s wing membranes were 
not suited to fast flight (36). 

In our analysis, we used the latest body mass estimates and 
drag coefficient estimates. For the drag coefficient, we used 
Palmer’s experimental results (34). For body mass, there is still 
variation even among recent studies (9, 15). Therefore, for Ptera- 
nodon, we employed the 2 estimates by Henderson (15) and Wit- 
ton (9). For Quetzalcoatlus, we only employed Witton’s estimate (9), 
since Witton and Habib (16) pointed out a problem with applying 
Henderson’s estimation method to giant azhdarchid pterosaurs. 











Models 


The dynamics of soaring animals are described using the equa- 
tions of motion (EOM). We first describe the EOM and parameters 
therein. We then describe the calculation procedure for dynamic 
soaring and thermal soaring, respectively. 





Aerodynamic forces and parameters 


We regard an animal as a point of mass. We employ the frame 
of reference (i, j, R) = (@east. north Cup) (44). Then, the dynamics of 
the animal’s 3D position X(t) = xi+ yj+zk and ground velocity 
V(t) = Vei + Vyj — Vsinkk are represented by the following EOMs: 


oe =L+D+mg (1) 
dt 
and 
dX 
errs =V. (2) 


When an animal is soaring, 3 forces—gravitation (mg), lift force 
L), and drag force (D)—act on it. Gravitation mg is a product of 
the constant of gravitation (g) and mass of the bird (m kg), and 
its direction is toward the ground. The direction of the lift force L 
is dorsal and perpendicular to the air velocity, Va (t) = V(t) — W(t), 
where W(t) represents the wind velocity. Drag force D is against 
the air velocity. For the analysis of dynamic soaring, we represent 
these EOMs in a different way by transforming the ground velocity 
to pitch y, yaw w, bank angle ¢, and airspeed V (= |Va|). The yaw 
angle w is the angle between i and the projection of Va to the 
ij-plane, and the pitch angle y is the angle between Va and the 
ij-plane and is positive nose down. We assumed that wind blow 
is along the y-axis and wind only depends on the altitude of the 
animal, i.e. W(t) = —W(z)j. Under this formulation the EOMs area 
represented by the following equations (44): 
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where W(z) represents the wind gradient. A specific form is pro- 
vided in the latter subsection. L represents the strength of the lift 
force, and D represents that of the drag force. The aerodynamic 


theory asserts that these values are 
1 
L=|L|= 5 PCL SwV" (9) 
and 


2 
D=|D|= 5 hCoproSwV? ae 5 CoparSBV" +e oe sv (10) 
Here, C, represents the lift coefficient, Sw represents the wing 
area, and p represents the air density and was set to p = 1.23 kg/m? 
(26). This is the International Standard Atmosphere value for sea 
level at 15°C expressed as 3 significant digits (1). We employed 
this value because, based on a geochemical model of the car- 
bon/sulfur cycle (79, 80), 20-23% of the air volume during the era 
of pterosaurs and giant extinct birds (approximately 70-6 Ma) was 
oxygen, which is almost the same as today (21%). We also tested 
the sensitivity of soaring performance to air density (See Air den- 
sity dependence of soaring performance). 
The drag is composed of 3 terms. The first term is the profile 
drag that stems from friction on the wing. Cppro is the profile drag 
coefficient. The Flight software developed for evaluating bird flight 
performance employs 0.014 as the default value for the Cppro (21). 
This value is often used for birds (1) and even for pterosaurs (16). 
However, theory predicts that the profile drag coefficient varies 
with the lift coefficient following 





Cppro = Cppro,min + Rpro (Cc. _ Cepek (11) 


and experimental results support this prediction (34, 81-83). 
Hence, for birds, we employed Cppro,min = 0.019, Rpro = 0.030, and 
C.ppromin = 0.77, estimated from a jackdaw gliding in a wind tun- 
nel (82). For pterosaurs, we obtained values from a wind tunnel 
experiment with reconstructed pterosaur wings (34). In the study 
(34), the reconstructed pterosaur wing span was approximately 
6 m; however, we assumed that the result would not vary sig- 
nificantly with the size of Quetzalcoatlus. We quantified the posi- 
tion of data points in Fig. 3(b) of (34) using WebPlotdgitzer and 
fitted (Eq. 11) to the data points (see Figure S3, Supplementary 
Material). In (34), results of 2 types of wing sections are shown: 
faired wings and unfaired wings. For faired wings, we obtained 
Cppromin = 0.050, Rpro = 0.11, and C, ppro.min = 1.0 and, for unfaired 
wings, Cppromin = 0.077, Rpro = 0.065, and Cy ppro.min = 0.89. Faired 
wings entail higher soaring efficiency compared to unfaired wings. 
For simplicity, we only show the results of faired wings in this 
study. 

The second term is the parasite drag stemming from friction 
on the body, where the Cpyar is the parasite drag coefficient, and 
Sp is the body frontal area. We used the following recently recom- 
mended formula: 


CpparSz = 0.01Sw, (12) 


on the practical basis that neither Cppar nor Sg is exactly known 
(45). 

The third term is the induced drag that stems from the lift 
force. Ra represents the aspect ratio (R, = b*/Sw, where b is the 
wingspan). The k is the induced drag factor; we set k to 1.1, as in 
previous studies (1, 16, 21). The lift coefficient has a maximum 
value; for birds, we set C;, to < 1.8 (21). As the aerodynamic prop- 
erties of pterosaurs can differ from those of birds, and the wind 
tunnel experiment indicated that C, could reach more than 2.0 
(34), we set the pterosaurs’ lift coefficient to < 2.0. 

The remaining parameters in the EOMs are body mass (m), 
wingspan (b), and wing area (Sw). For these morphological param- 
eters of extant birds, we used values reported in previous studies 
(shown in Table 1). 
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The EOMs include variables that soaring animals can control, 
that is, bank angle ¢(t) and lift coefficient C;. Although these vari- 
ables are time-dependent, for simplicity, we assumed that the an- 
imals keep their lift coefficients at a constant value. Hence, using 
a time series for bank angle, a constant value of C,, and values 
of parameters, the dynamics of the soaring animals were deter- 
mined with EOM. 





Quantification of the dynamic soaring 
performance and the minimum required wind 
speed 

Wind gradient models 


We explored 2 types of wind gradients. The first was the logarith- 
mic model represented as 
[Logarithmic wind gradient model] 


Wio 


Mes!) = nLOfhnin) 


In(Z/hmin)- (13) 
This function is defined at z > hin. We set Nin to 0.03 [m], fol- 
lowing a previous study (40). W1o is the wind speed at height z = 
10 m. This model is deemed to be a good model of the average 
wind field in the first 20 m above the sea surface, assuming a flat 
sea surface, and has been a popular approach in dynamic soaring 
modeling. However, recent studies argued that the real sea sur- 
face is not flat, and wind separations in ocean waves may occur 
more often than expected (60). To describe wind-separation-like 
wind profiles, a sigmoidal model has been proposed (44, 84). We 
also employed the sigmoidal wind model with a minor change, 
represented as 

[Sigmoidal wind gradient model] 
Wax 


Wsigmoid (z) = hy * 
1+e7 — 








(14) 





The hy, determines the height of wind separation, as shown in 
Fig. 3. In this study, we set hy to 1, 3, and 5. The 4 is the thickness 
parameter. The wind speed changes with height (|z — h,| < 36 m). 
In a previous study, the wind shear thickness was speculated as 
approximately 1.5-7 m. Here, we set 6 to 3/6 with a steep wind 
change, and 7/6 with a gentler change (Fig. 3). 


Formulation to numerical optimization 

The numerical computation of dynamic soaring performance and 
minimum wind speed boiled down to the restricted optimization 
problem (85). That is, a mathematical problem to find the values 
of (i) a certain variable Y that maximizes (ii) an objective func- 
tion f(Y), satisfying (i1i) equalities h(Y) = 0 and (iv) inequalities g(Y) 
< 0. In the following, we describe the variables, object functions, 
equalities, and inequalities for dynamic soaring. 


Variables 


The dynamics of dynamic soaring animals are described by the 
D position (x(t), y(t), z(t)), pitch angle y(t), yaw angle w(t), air- 
peed V(t), bank angle g(t), lift coefficient C_, and the period of 
1 dynamic soaring cycle tr. Among these variables, 3D position, 
pitch, yaw, bank, and airspeed are functions of time t (0<t <r). 
Optimization problems that include functions as variables are dif- 
ficult to be directly solved. Therefore, we employed a collocation 
approach (43, 44). The collocation approach discretizes the vari- 
ables in time, such as X(t) (0 < t < r) to variables X; = X((i-1)N/r) 
(i= 1, N), and converts the EOM to the equalities between those 
discretized variables. Hereafter, we use X1:" = {Xi, X2,..., Xy}. In 
this study, we set the number of discretization points to N = 51in 
order to perform computations with reasonable accuracy within 


nA W 











a reasonable amount of time. Accordingly, the variables of this op- 
timization problem are position X1:n, yin, Z1.n, pitch angley:.n, yaw 
angle Win, airspeed Vin, bank angle ¢y.n, lift coefficient C,, and a 
period of 1 soaring cycle t. In addition, when computing the min- 
imum wind speed required for sustainable dynamic soaring, Wio 
log model) or Wmax (sigmoid model) were also treated as vari- 
ables. Hence, the total number of variables were 7 x 514 2 (+1 
when computing the minimum wind speed]) = 359 (or 360). 





Object function 


First, we computed (1) the minimum wind speed required for sus- 
tainable dynamic soaring for each wind gradient model. As the 
objective function to minimize, we set Wi for the logarithmic 
model and Wmax for the sigmoidal model. Then, we computed (2) 
the maximum travel speed averaged over 1 dynamic soaring cy- 
cle by maximizing the object function to ,/x%, + y%/t. Finally, we 
computed (3) the maximum upwind speed averaged over 1 dy- 
namic soaring cycle by maximizing the object function yy/r. With 
respect to the maximum travel speed and maximum windward 
speed, we computed these values for different wind speeds, that is, 
from the minimum required wind speed of the species to the high- 
est minimum required wind speed among the examined species 
(i.e. Quetzalcoatlus) + 2 m/s. In this wind speed range, the maxi- 
mum travel speed reached an unrealistically high value and/or 
the optimization calculation did not converge for some species. 
Thus, we stopped the computation of the maximum travel speed 
at the wind speed where the maximum travel speed exceeded 
40 m/s (144 km/h). 








Equalities 





The first equalities to be fulfilled for dynamic soaring animals are 
given in Eqs. 3-8. The collocation approach converts the EOM into 
the equalities between the variables listed in the above section. As 
the number of original EOM was 6 and the number of discretiza- 
tion was 51, the EOM were converted into 6 x 51 = 306 equalities 
(see (43, 44) for the specific representations of these equalities). 

The second type of equalities to be fulfilled were periodic 
boundary conditions of dynamic soaring: at the beginning and end 
of 1 dynamic soaring cycle, the state of the animal (..e. pitch, yaw, 
airspeed, bank, and height) is the same, represented as 





Z1 = Zn, ¥1 = yn. V1 = Wn, 61 = On, V1 = Vy. (15) 





Inequalities 

First, we assumed that there was a maximum limit of physical 
load on the animal. This is because dynamic soaring entails dy- 
namic maneuvering, which results in a corresponding accelera- 
tion. We employed the approach of a previous study (40) that re- 
stricted the load factor (L/mg) to less than 3, 


L 
— ie 3: (16) 

mg 
The second inequality was an important modification of the 
previous models. The height of the animal’s wingtip above the sea 


surface (Zwing) was calculated and represented as 
b. 
Zwing = Z— 5/sin | cosy > 0. (17) 


Previous studies discarded the existence of the sea surface (44) 
or restricted birds to only flying higher than a given height (1.5 m) 
from the sea surface (40). However, the height an animal can fly 
depends on the wing length and the bank angle (e.g. with a shorter 
wing length and a lower bank angle, an animal can fly at a lower 
height). When dynamic soaring birds fly, they adjust their wingtips 
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close to, but avoid touching, the sea surface (Fig. 3D). Dynamic 
soaring animals can exploit more flight energy when they pass 
through stronger wind speed gradients. As the wind speed differ- 
ence is strong close to the sea surface, how close to the sea sur- 
face an animal can fly is crucial for dynamic soaring animals. Ac- 
cordingly, long wings may restrict the minimum height at which 
the animal can fly and disturb efficient dynamic soaring. Hence, 
considering the effect of wings is crucial for evaluating dynamic 
soaring performances. 

Third, we assumed that the height of the animal was higher 
than 0.5 m, that is 


z>0.5. (18) 


The optimization problem described here is a restricted nonlin- 
ear optimization problem. We used the SQP method to solve the 
problem with the “fmincon” function in MATLAB® Ver R2019a. 


Quantification of the thermal soaring 
performance and the required upward wind 
speed 

For the computation of glide polars and circling envelopes, we fol- 
lowed a similar procedure to the Flight software, as described in 
(21), but modified some parameters and settings. Below, we out- 
line the procedure and parameters employed in this study. 

First, to compute the glide polars and circling envelopes, we 
needed to set rules regarding how gliding animals adjust their 
wingspan with respect to their airspeed. To achieve this, we used 
2 wingspan adjustments: linear wingspan reduction and fixed 
wingspan (21). Although a fixed wingspan resulted in a better 
performance than linear wingspan reduction, we found no sub- 
stantial differences in the results (see Fig. 5; Figures S4-S6 and 
Table S1, Supplementary Material). Hence, for birds, Fig. 5 shows 
the results of adopting linear wingspan reduction, which is Flight’s 
default setting based on measurement of a jackdaw gliding in 
a wind tunnel (86). The default settings of Flight assume that 
wingspan, wing area, and thus aspect ratio linearly decreases with 
factor B = (Bstop- V/Vs)/(Bstop- 1). Therefore, we replaced b, Sw, and 
Ra with Bb, BSw, and BRa, respectively. In this equation, Vs is the 
stall speed, the airspeed of the animal at the highest lift coeffi- 
cient (i.e. Vs = /(2mg)/(eSwCimax)); and Bstop is a constant that de- 
termines the degree of wing reduction. We set Bestop to 5, the de- 
fault value in Flight. For pterosaurs, we show in Fig. 5 the results 
of a fixed wingspan. This is because, to the best of our knowledge, 
no study has yet shown a clear mechanism, whereby a membrane 
wing can be reduced in span and area without slacking the mem- 
brane. For reference, we also show results of the linear wingspan 
reduction, where Bestop was set to 6, in Figures S4-S6 (Supplemen- 
tary Material) and Table S1 (Supplementary Material). This is a 
value employed in a previous study (16), although no quantitative 
evidence for this value was provided therein. 

Then, the glide polars were derived from the EOM, setting bank 
angle to 0, assuming that the pitch angle was small enough (y « 
1) and considering the gliding animal was at kinematical equilib- 
rium (mg = VL? + D? ~ L, siny = D/L ~ D/mg). Sinking speed was 
represented as a function of airspeed V (21), 


wee ia if 1 2k ) 1 
—— 2 (=) Perea (2) (sos Vv" a) 








Note that we used Eq. 12 (Cppar Sg = 0.01Sw; in this equation, Sw 
is not replaced with BSy) to derive the above equation. Eqs. 11, 19, 
and the equilibrium relation, C, = 2mg/(efBSwV7), result in the 


following equation that gives a glide polar: 


S 2m 2 
Vsink = 5 (=) f conn + Rpro (> _ C..dpein 


3, (mg 2k 1 
soon}? (88) (ss) a 


The horizontal speed is ,/V — Vé,,. Thus, the maximum glide 





ratio is the maximum value of ,/V?2 — Voink /Voink- 

Next, we examined the circling envelope. We considered a sit- 
uation where an animal is gliding on a circular track of radius r 
with a bank angle of @ and the lift coefficient C*. The circle radius 
and the sinking speed in that situation (Vsink circle) 18 given by the 
airspeed V(C;) and the sinking speed Vsink(C7) when the animal 
is gliding along a straight track with a lift coefficient C7 as below 
(87). 








vcr 
= (i) (21) 
gsing 
and 
Vsink (Cf) 
Vsink,Circle = (cosa)? ay? , (22) 
where V (C7) = ,/2mg/(pBSwCF), and Vsink (C7) is obtained by substi- 


tuting V(C}) into Eq. 20. 

There are 2 types of what we refer to as a "circling envelope’ (21, 
22, 88, 89). The circling envelope, proposed by Pennycuick (22), is 
an upper performance limit that an animal can achieve in the as- 
cending phase of thermal soaring. The circling envelope was de- 
fined as a function that gives, for a circle radius r, a minimum 
sinking speed that an animal can achieve when gliding in a circle. 
Specifically, this is expressed as 





MING <Cimnax.0<$<¥ Vsink,circle(V (Cz) , @) with the constraint 


vcr 
~ gsing 





(23) 


At the limit where r > oo (.e. ¢ > 0°), this represents straight 
gliding flight with minimum sinking speed. The Cf on the circling 
envelope varies with radius, and in the small radius region, it 
reaches to the maximum lift coefficient. This means that the an- 
imal is flying at a stall speed. Note that r should be larger than 
V(Cimax)*/g, which is the limit of CF = Crmax and @ > 90°. This min- 
imum value is referred to as the limiting radius, i.e. the animal 
cannot achieve a circle radius smaller than this. 

Furthermore, what Pennycuick newly called a circling envelope 
was a plot of Eqs. 21 and 22 as parametric equations with ¢ as the 
mediating variable, assuming C} to be constant (21). As for Cy, the 
lift coefficient minimizing the sinking speed during a straight glide 
was employed (21) (although there is a study that employed the 
maximum lift coefficient (90)). In this definition, at the limit where 
¢ — 0°, as in the former circling envelope, the animal is gliding in 
a straight line with a minimum sinking speed, but the limiting ra- 
dius for @ > 90° is greater than the former. This is because the 
latter circling envelope assumes that the animal does not change 
its lift coefficient regardless of the circle radius. In other words, 
the former circling envelope shows an upper limit of the ani- 
mal’s capabilities, whereas the latter implicitly assumes that the 
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animal is circling with some margin of safety that avoids the risk 
of stalling. 

In recent years, the latter has often been used to determine the 
circling envelope (21, 51, 90). However, as the wings of pterosaurs 
are expected to have been able to achieve a higher maximum 
lift coefficient than that of birds (34), pterosaurs might have had 
a greater advantage than birds in increasing the lift coefficient 
when circling with a narrow radius. Therefore, to account for these 
differences, we evaluated both circling envelope definitions in this 
study. The 2 circling envelopes differed in the small circle radius, 
but there was no significant change in the relative performance 
relationship between the species. See Fig. 5 and Figure S5 (Supple- 
mentary Material) for results of the former and Figure S6 (Supple- 
mentary Material) for the latter results. The limiting radius, min- 
imum sinking speed in straight gliding, maximum glide ratio, and 
horizontal speed at maximum glide ratio are shown in Table S1 
(Supplementary Material) for the 2 wing reduction methods and 
2 circling envelope definition. For ASK 14, the measured glide po- 
lar (Vsink = 42 + 5.2 x 10-°V3) and the circling envelope of the for- 
mer definition setting maximum lift coefficient to 1.3, reported in 
a previous study (22), are presented. 














Air density dependence of soaring performance 


To test the sensitivity to the air density, we also evaluated the soar- 
ing performance in high air density environments that pterosaurs 
and giant birds would have never experienced, 1.e. at an air den- 
sity of 1.2 times the current level (1.48 kg/m?). This air density 
corresponds to the value that would be obtained if the oxygen 
concentration increased to 30% and the amount of other compo- 
nents were equal to those in the present day (91). This air den- 
sity is also expected to be higher than that during the eras of 
pterosaurs and giant extinct birds, as the concentration of oxy- 
gen reached 30%, approximately 300 Ma (Permian period), and 
had already decreased to a value comparable to modern day air 
in the eras of giant pterosaurs and giant extinct birds (70-6 Ma) 
(80). 

We evaluated the dynamic soaring performances of P. sandersi 
and Pteranodon under this higher air density (see Figure S7, Sup- 
plementary Material) with wind gradient shown in Fig. 3. Al- 
though the higher air density resulted in a slightly lower mini- 
mum wind speed required for dynamic soaring, the qualitative 
results were similar to those obtained with the modern day air 
density (see Figure S7, Supplementary Material). We also evalu- 
ated the thermal soaring performance of Quetzalcoatlus under the 
higher air density (see Figure S8, Supplementary Material). De- 
spite the high air density, the thermal soaring capability of Quetzal- 
coatlus remained low. Overall, our results were robust to air density 
uncertainty. 
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